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Introduction. 


Starting from the first paper of E.Witten [2] of 1981, the method of Supersymmetric 
Quantum Mechanics (SUSY QM) became a new effective tool [3], [1] for investigation of 
different problems of conventional Quantum Mechanics. From a mathematical point of 
view, this method - for one-dimensional space - is a reformulation of so called Darboux 
transformation of Sturm-Liuvillc equation well known to mathematicians since more than 
hundred years [5]. SUSY QM approach was helpful for investigation of different problems 
in Quantum Mechanics and a lot of generalizations were proposed. In particular, the multi¬ 
dimensional generalization i, 0. IS of SUSY QM which led to many new results, mostly 
for two-dimensional case 0. ra has to be mentioned. This multi-dimensional approach can 
be considered also as a generalization of old Darboux transformations. 

Among new ideas provided by SUSY QM method, one - the shape invariance - is of 
special significance 0, n, mi, mi. It allows to solve some quantum models exactly or 
quasi-exactly. The list of known exactly solvable one-dimensional models is very restricted 
(approximately, about ten potentials are in this list). Supersymmetric transformations pro¬ 
vide the partnership between pairs of potentials and therefore lead to new solvable potentials 
ra, though sometimes of rather complicated analytical form. The method of shape invari¬ 
ance gives a very elegant algebraic algorithm to construct exactly solvable models without 
direct solution of differential equations. All previously known [Qj one-dimensional exactly 
solvable systems were shown to be shape invariant ones pg. Also, this approach provides 
analytical solution of several nontrivial two-dimensional problems [10], PH. 

Recently, the general investigation of additive form of shape invariance for superpotentials 
without explicit dependence on parameters was performed by J.Bougie et al p, pn. it 
was demonstrated that no additional shape invariant models in one-dimensional case can 
be constructed. This result was obtained in the framework of standard SUSY QM with 
supercharges of first order in derivatives. Meanwhile, it is known [18], [19], [20] that such 
standard SUSY QM does not exhaust all opportunities to fulfill the generalized SUSY algebra 
for Superhamiltonian H and supercharges Q ± . Higher order supercharges are also possible, 
and for example supercharges of second order in momentum lead to some new results [20] . 
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In the present paper just a potential which fulfills the second order supersymmetry will 
be produced by the condition of shape invariance of second order. This potential is of very 


compact form, it depends on arbitrary parameter but has a strong singularity of g/x 2 kind 


with g G (—1/4,0). This last property forces to construct a suitable self-adjoint extension 
of the Hamiltonian H, i.e. the suitable class of functions where H acts. The paper is 
organized as follows. Section 2 contains a brief summary of one-dimensional SUSY QM and 
shape invariance. A shape invariant potential is built in Section 3 by means of second order 
supercharges. The direct analytical solution of the Schrodinger equation is given in Section 
4. The problem of a suitable class of functions belonging to the domain of the self-adjoint 
extension of H is discussed in Section 5 in terms of connection conditions at the origin. 
Section 6 includes derivation of the Spectrum Generating equation and its analysis, and 
finally, some conclusions are given in Section 7. 

2. SUSY QM and Shape Invariance. 

The main ingredients of one-dimensional SUSY Quantum Mechanics and of shape in¬ 
variance will be briefly presented in this section (see details in Refs.[3]). An arbitrary one¬ 
dimensional Hermitian Schrodinger Hamiltonian with potential V(x) and mass m = 1/2 can 
be factorized: 


H = Q + Q~ = -h 2 d 2 + V(x)\ d=4~ 

dx 


( 1 ) 


by means of the first order differential operators Q ± ; Q + = (Q )'. Directly by construction, 
the superpartner Hamiltonian 


H = Q~Q+ = -h 2 d 2 + V{x) 


( 2 ) 


is intertwined with initial H : 


HQ + = Q + H ; Q~H = HQ~. 


(3) 


Here, the first order supercharge operators Q ± are expressed in terms of the superpotential 
function W(x) : 
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and the partner potentials are: 

V(x) = W 2 (x) - hW'(x)] V(x) = W 2 (x) + hW'(x). (5) 

As a consequence of (j3jl . almost all wave functions of //. H are interrelated: 

^n(x) = Q + ^ n (x)-, ^n(x) = Q~^ n (x), (6) 

the only exclusions are possible normalizable zero modes of operators Q ± , but each of oper¬ 
ators admits not more than one such zero mode. 

In some cases the partner Hamiltonians H , H have more close relationship. Let us sup¬ 
pose i, n, mi, im [16] that both Hamiltonians ([I]), potentials (J5]), superpotentials depend 
not only on coordinate x, but also on some parameter a, something like a coupling constant: 
H = H{x\a ); H = H(x]a), etc. Let us suppose also that the potentials V(x]a), V(x,a ) 
obey the so called shape invariance property: 

V(x; a + h) = V(x; a). (7) 

Mathematically, the Planck constant h above can be considered above as a constant with arbi¬ 
trary numerical value. Practically, in the Schrodinger equation context, the shape invariance 
of known solvable potentials is realized just by h shift of parameter (see details in [1], pTJ). 
In essence, the relation ta means that the dependence of superpartners H(x]a), H(x;a ) 
on coordinate x is qualitatively close, but quantitatively slightly different. The distinction 
is determined by dependence of potentials V, V on slightly different values of parameter, 
i.e. ’’the shape” of potentials is the same. The potentials V(x;a), which obeys the relation 
o. are called shape invariant potentials. The list of such potentials is very short, and it 
practically coincides with the list of known solvable potentials. Usually, exact solvability 
of the Schrodinger equation with specific potential is established by exact solution of the 
corresponding Sturm-Liuville equation providing the analytic expressions for normalizable 
eigenfunctions and eigenvalues - the wave functions and energies of quantum system. In con¬ 
trast to this, shape invariance (JTJ) together with intertwining relations (j3j) provide different 
- purely algebraic - algorithm to find the spectrum and wave functions of the model. 

The algebraic shape invariant algorithm gives the following prescription. If one chooses 
the ground state energy of H(x, a ) to be E 0 (a), the energy of n— th excited state of H(x, a) 
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is: 

E n (a) = E 0 (a + nh), 

and its wave function is obtained from the ground state wave function as 

a) = Q + {x , a)Q + (x, a + h)...Q + {x, a + (n — l)h)tyo(x, a + nh). (8) 

The ground state function T 0 (a;, a) is supposed to be known for arbitrary values of parameter 
a. The shape invariance relation (J7J) and intertwining relations (E]) provide, for example: 

H(x, a)Q + (x, a)T 0 (x, a + K) = Q + {x, a)H(x, a)T 0 (x, a + K) = 

= Q + (x, a)H(x , a + h)^ 0 (x, a + h) = E 0 (a + h)Q + (x , a)T 0 (x, a + h). 

This chain of relations, obtained purely algebraically, shows that indeed, is an 

eigenfunction of H(x,a ) with energy E 0 (a + K). Analogously, the expression OSJ) can be 
checked for arbitrary natural n. 

In papers [11] the most general dependence of the superpotential in (J4]) on parameter a 
was studied. It was proved that all conventional exactly solvable potentials, such as harmonic 
oscillator, Coulomb, Morse, Poschl-Teller, Scarf etc. can be reproduced by means of the 
shape invariance approach with W = W(x, a). Meanwhile, some recently constructed exactly 
solvable potentials ra were missed in this algorithm. In paper [1], a more general class of 
supercharges was considered, where the superpotential depends not only on parameter a, but 
also explicitly on the constant h. Expanding the superpotential in powers of h , potentials 
of PE] were also obtained as a particular case. In the next Section, the approach with this 
more general dependence of supercharges will be used, but in the context of higher order 
SUSY QM. 

3. Construction of New Shape Invariant Potential. 


The conventional SUSY QM admits the so called Higher Order generalization [18], [19]. 
ra where the operators Q ± are differential operators of order IV > 1. It is known that 
arbitrary higher order transformation is equivalent to a series of first and second order 
elements. In turn, the latter can be factorized onto two first order supercharges which 
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intertwine the initial Hamiltonian with the intermediate Schrodinger-like Hamiltonian [19]. 
The intermediate system has real or complex potential depending on the sign of the constant 
parameter d (see (J9]) - (Hill below), and in this sense second order supercharge is [19j reducible 
or irreducible, correspondingly. 

The case of supercharges of second order was studied completely na, a, and just this 
case will be used in the present paper in the shape invariance context. The most general 
form of second order operators Q ± is: 

Q + (x, a) = h 2 d 2 — 2 hf(x, a)d + b(x, a); d = d/dx (9) 

It is known that the intertwining relations (j3j) with supercharges (J9]) are fulfilled iff 

n*,a) = -2 Hf + f + ^-^)-XL +1 , (10, 

where f(x) is an arbitrary function of x and function b(x) is expressed in terms of f(x) : 

b(x) = -hf + f - + (ftT) +_1 (11 ) 

Both d and 7 are arbitrary constants depending explicitly only on parameters a and h, in 
contrast to f(x), which depends only on coordinate x. This form of supercharges Q ± is not 
the most general, but the choice of a—dependent functions f(x) would lead to much more 
complicated problems. 

Thus, the case with 

OO OO 

f(x, a) = f(x), d(a, h) = ^ h n d n (a), 7(0, h) = ^ h n ^ n (a) 

71=0 71=0 

will be considered. Resummation of series gives: 

00 n—1 \ 00 

d(a + h,h)~ d(a, h) = y rr V y = V h n D n {a), 

in — k)\ ' 

n= 1 k =0 ' ' n= 1 

00 n—1 ( n—k )/ \ oo 

7 (a + ft, ft) - 7(o. fi) = E fi " E fa _ mi = E h " A M' 

n=l fc=0 ' n=l 

Then, the shape invariance condition ([7]) leads to relation: 

4hf( x ) + S ”- 1 p ( fa ,l(a) = E 

' n=l 
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It is obvious now that: 


A f'( x ) + 4 /a (a;) “ Al ^ 

D n (a ) = A n (a) = 0, n > 1. 


( 12 ) 

(13) 


In turn, taking derivative of (1T2|) over parameter a and separating variables, one obtains for 
nontrivial f(x) : 


DA a) = const = 


16c 3 


Ax (a) = const = a 


(14) 


where the definition of the constant for Di in the r.h.s. will be clear below (see Eq. ffTT]) ). 
Eqs. (IT3l) . (fT4l) give: 

d n (a) = const = d n , 7 n (a) = const = j n , n > 1 - 
By definition, D\{a) = d' 0 (a), Ai(a) = 70 (a), and the dependence has the form: 

-3 _ 22 _ 4 c 3 

i+K{fi), 7 ( 0 , h) = aa-\- 

n =1 


d(a, h) = 


16c 3 


-a- 


_^ 4c° _> 

-^2h n d n =-—a + n(h), 7 ( 0 , ft) = cra + ^ h n ^ n = cra + <5(ft), (15) 


n=l 


where the constants c, k, cr, 5 do not depend on the parameter a. 

Thus, the potential (ITU with constants d, 7 from (fT5]i is shape invariant if the function 
f(x) satisfies equation: 


/'(*) 


= a/A. 


(16) 


3 P(x) 

The case of c = 0 in (1I6|) corresponds to the well known potential - singular oscillator 
potential - and it will not be considered here. For all nonzero values of constant a, equation 
(TT6l) can not be solved explicitly. This is a reason to choose the a = 0, for which one obtains 
the monomial solution: 

f(x) = or 1/3 , (17) 

and, after the convenient choice ft = 1 (it can be fixed independently from the earlier choice 
2 m = 1 ), the superpartner potentials are as follows: 




( 18 ) 










It is easy to check that by rescaling the coordinate x —» c~ 3 ^ 4 x in both kinetic and potential 
terms of the Hamiltonian one can achieve c = 1. From now on this value of c will be used, 
and thus the model depends only on one parameter a. Potential (1T8]l is even function of x, 
its form is shown on Figs. 1-3 for different values of parameter: a = 1, a = 2 and a = 3. 



-6 -4 -2 2 4 6 


Fig.l Plot of potential (Tl8|) for a — 1. 



Fig.2 Plot of potential (H8]l for a = 2. 
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Fig. 3 Plot of potential (1T8|) for a = 3. 


It is obvious that, the potential V(x) is shape invariant in parameter a, as it was expected 
by construction: 


V(x; a + 1 ) = V(x] a); 7 (a) = 0 . 


(19) 


As a remark, from the very beginning one could avoid the dependence of d and 7 on h , and 


to check directly that the dependence (TT5]i of these constants only on a provides both the 
intertwining relations and shape invariance. 

4. Direct Solution of the Schrodinger equation. 

Our task is to investigate the potential V{x) in (jT8|l for positive values of constant 
a > 0. This potential has a strong but physically acceptable singularity ~ \fx 2 at the 
origin: its coefficient is larger than —1/4 and prevents the ’’fall to the center” [21]. In this 
Section we will demonstrate that the Schrodinger equation with this potential can be solved 
analytically in terms of confluent hypergeometric functions [22]. Due to the singularity at 
x — 0 , this solution is not sufficient to provide the physical wave functions of the problem. 
Indeed, while the growing behaviour of potential (TT8]i at infinity corresponds to the unique 
decreasing (L 2 —integrable) solution, the situation at the origin is different. Namely, both 
independent solutions at the origin are decreasing, and this fact does not allow to determine 
the discrete spectrum of the model. According to standard terminology [23], the behaviour 
at x —> i00 is of the limit-point type, while the behaviour at x —> 0 is of the limit-circle 
one. How to deal with this problem will be studied in the next Section of the paper. 

Let us look for the solution of the Schrodinger equation with potential (fT 8 j) : 


H{x)V{x) = EV{x)] H(x) = -d 2 + V(x) 


( 20 ) 


in the following form: 



( 21 ) 
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where according to (ITT]) f{x) = x 1 / 3 . Choosing A (a) = — the equation (l20j) can be 


rewritten as: 


w , S'{x) ~ 2 P(x) + A(a) 


d z + 


f(x) 


d + E — A (a) ) il>(x) = 0. 


In the following, the variable 


C = V3(f(x) - -) 


will be used instead of x, and Eq. (122]) takes the form: 


d 2 d \/3 . „ ,. ^ 3A . ^ ,, . . 

d? “ C dc + X 1[E - A)c + T {e - A) hK) = 


Separating one more multiplier from "0(C) : 


- a/3 


V’(C) = 2/(0 exp (—(£-A)C), 


one obtains new form of Eq. 


i? + ( ir (B " A) “ c) ^ + 0 £2 " A2) ) y(c) “ °- 

Finally, the new transformation of variable: 

( = \(C-^E-X)) 2 = 3 -(2f( x )-E) 2 . 

reduces Eq. fEZBD to the equation for the confluent hypergeometric function: 

^ + q - 4 (a2 - £2) - 

It has two independent solutions: 

</iR) = e(«, l;f); »«) = C I/2 e(« +1 |;0, 


( 22 ) 

(23) 

(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


with F - confluent hypergeometric functions H- The general solution of initial equation 
(120|) is a linear combination (iV lj2 —arbitrary constants): 


*(i,S) = Af,4 (1) (i,B) + Af 2 t (2) (i,B), 


(30) 


with 


= T 0 (x, a) exp (^-(E - A )Qyi(£{x)) ~ ^o(x, a)yi(£(x)) exp - A)/ 2 (x) ). (31) 
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The asymptotic behaviour of the confluent hypergeometric function for positive £ —> oo 


is: 


F(a, c; 0 ~ 

T(a) 


(32) 


where T is the standard Gamma function. Thus, taking into account the x —> ±cx) behaviour 
of all multipliers in (1ST]) , the asymptotics of T is: 

p4 o f4 T-1/Q /q\ q/4 q/4 




3/ 4 3/ 1 1/2 

ex P (—)(—) ' 


• £C1^ + JV2 _ 2v / 2r(3/2) 


r(a) 


V3r(a + 1 / 2 ) J ’ 


(33) 


where some constants were absorbed additionally in normalization constants Ni,N 2 . The 
condition of vanishing of the exponentially growing terms in asymptotics provides relation 
between constants , N 2 : 


Ni _ T(3/2) T(q) 2^2 = _ ygT(a) 

iV 2 ~~ T(l/2) T(a + 1/2) ^3 “ 7 _ 73T(a + 1/2) ’ 

and one of them is still arbitrary (up to normalization of wave function (j30]l ). 


5. Connection conditions at the limit-circle singular point. 


Thus, in the previous Section, the general L 2 — integrable solution (130)1 . f)3T]l of the 
Schrodinger equation for arbitrary value of energy E was obtained up to an overall nor¬ 
malization constant. This knowledge does not provide the spectrum of the model: the 
energy E does not appear to be constrained, while from the form of potential ()T8j) it is clear 
that the spectrum should be discrete, bounded from below and not bounded from above. 
The problem is well known since a long time [23] in the context of different physical models, 
and its origin is due to the properties of the Hamiltonian with potential dT8|) around the 
singular point x —» 0, which are reflected in the behaviour of solution (1301) in this region. 

The point x — 0, where both independent solutions are normalizable, is the so called 
limit-circle point. As for Hamiltonian of the model, it is not self-adjoint, and one has to find 
some self-adjoint extensions of H, which are acceptable from physical and/or mathematical 
point of view. The task of such sort was considered in many papers, all of them describe 
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the algorithms to find the self-adjoint extensions. It was proved that the entire variety of 
extensions forms the two-parameter set described by elements of unitary group 17(2). We 
will follow one of such algorithms - due to papers IS! - formulated in terms of connection 
conditions in the limit-circle singular point (see also the preceding papers [26] . 

Below the scheme of construction of the entire variety of self-adjoint extensions of given 
Hamiltonian with singularity of limit-circle type at the origin will be briefly presented. The 
first step of the procedure is to construct a pair </?W(x) of independent auxiliary solutions 
(reference states) of the Schrodinger equation (1201) with some (arbitrary) real energy value e. 
The normalizability of reference states at infinity is not assumed, and their behaviour at the 
origin is determined by the singularity of potential - in our case of potential (fl8l) independent 
solutions have the following behaviour at the origin: ~ x 1 ^ and ~ x 5 / 6 . 

Let us suppose that the auxiliary solutions of (1201) are known, and that their 

Wronskian is normalized: 


w[(pW,<pV>] = tpWtpW - ^V 2)/ = i. 


Let T(a;) and T (x) to be arbitrary decreasing at infinity functions for which differential 
operator is well defined. The condition of symmetric property of H reads: 


dx['^*(x)H(x)^ f (x)] = / dx[(H(x)^/(x)y^/(x)]. 

J x 


(35) 


Integrating (l35j) by parts and taking into account that both T and T are exponentially 
decreasing at infinity, one obtains ’’the boundary terms” only from the left and right vicinities 
of x = 0 : 

0 = J *:[f*(i)ff(i)*(a:)] - J dx[(H(x)'i(x)Y'i(x)] = (w^t*, U+o - W\¥, <t]_ 0 ) ■ 

(36) 

The r.h.s. of (j36j) can be expressed using the auxiliary solutions </?7) v i a the chain of the 
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following transformations of Wronskians: 




vj/* vj/*' 


vp* 



T T' 


T 

T' 



<pW 

-<pW 


cpW 


VtpW - tv (1) ^ (2)/ - ^V (2) 

= Vk[^*,(^ (1) ]Vk[^,^ (2) ] - W[ty*,ip {2) ]W[V,<p {1) ]. 


This transformation allows to rewrite the r.h.s. of (136|) as: 




(37) 


where the two-component columns are defined as: 

nW =/ w\*,vmu y n(2) / ^[*.^>1+0 \ 

( W'[*,P (1) ]-o J ’ \ -W[' I,v> |2) ]-o ) ’ 

and analogously, are defined by replacing T —> T. For T and T in a self-adjoint domain 
the expression (j37j) has to be zero. This condition can be rewritten in terms of 

= h2 ( ' 1 ^ ± (38) 

with an arbitrary nonzero constant u. For coinciding T = T, the vanishing of (1371) is equiv¬ 
alent to the relation: 

= UeU{ 2 ). (39) 

For the general case of different T and T, but with the same matrix [/, the difference (1371) 
actually vanishes, and thus the symmetric self-adjoint operator H is completely character¬ 
ized by matrix U, which may be therefore called characteristic matrix of the self-adjoint 
Hamiltonian. 

The scheme described above can be directly applied to our model with potential (fT 8 j) , 
starting from the convenient choice of the reference modes = 1,2: 

= — ^$^(x, e = A) = — ^\l/ 0 (x, A) forx> 0; ^^(— x) = — ^^(x); (40) 

< P {2) = exp (—4a)<F ( 2 ) (o:, e = A) = exp (-4a) (2f 2 (x) - X)^ 0 (x, \)F(l/2, 3/2; £); 










where constant multipliers were chosen so to provide the value of Wronskian: W[^ l \ = 

1. It will be clear below that the freedom in this choice does not affect the final results. 
Indeed, for diagonal U equations for spectrum of the next Section does not depend on these 
constants, and for nondiagonal U the change of constants can be compensated by change of 
ui in ([38]) . 

The next step is to choose the specific matrix U G 17(2), which defines the symmetric 
operator H. For example, the class of operators with diagonal U can be selected. Namely, 
physically it corresponds to the vanishing probability flow through the point x — 0 , since for 
such U, the boundary conditions at x —> +0 are not mixed with those for x —> —0 : 

j(x = 0) = — = 

2 inn 

= —(W[V*, </> (1) ]IT[T, </> (2) ] - </> (2) ]IT[T, </> (1) A = 0. 

2 im \ ) 

Below the simplest ansatzes U = ±7 will be considered. For the choice U = —I, the 
first two-component column has to vanish f7') = 0. Since both solutions $7) _£;) are even 
functions of x and </7*) are also of definite parity, the components in Q 1 - 1 ) differ only by sign 
from each other. Thus, the Hamiltonian H is self-adjoint in the domain of functions T 
satisfying the only - on the right side - condition: 

W[f(x,E),<pW(x)} +o = 0, (42) 

which realizes the connection condition for 4/(a;, E) in the limit-circle singular point. This 
equation depends on energy E, and it will determine the allowed values E n , i.e. the discrete 
spectrum of the model. The wave functions (1301) . specified by fl3Tf and reference mode (1401) 
must be substituted into (1421) : 

Nt • lT[<F ( 1 ) (a:,E n ),^ ( 1 ) (a ;)] +0 + N 2 ■ W [<h (2) (x, E n ), </? (1) (®)] +0 = 0. (43) 

Due to x —y ±oo asymptotic, one of the constants Ni,N 2 is expressed in terms of second 
(see (JMD), which is still arbitrary up to normalization. 
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6. Derivation of the Spectrum Generating Equation. 


Practically, to find the Spectrum Generating Equation (SGE) for U = — I one has to 
rewrite Eq. (14 3 P by substituting functions 'hi*' and y/ 1 ) from d3T]) and (HOI) . The expression 
(/jd) allows to obtain the x —> +0 limit of Wronskians in (HOI) explicitly as linear combinations 
of Confluent Hypergeometric Functions: 


W[& 1 \x,E),<p( 1 \x)] +0 = -- 


—AaEF(a + 1, 3/2; 3E 2 /8) 


(E — X)F(a, 1/2; —- 


W{&V(x, E), f B(4, = - 3E{E S X) F(a + 1/2,3/2; 

, , , 3E 2 . 3E 2 (2a +1) , , , 3 E 2 s 

~F(a + 1/2, 3/2; —)- K — - ’-F{a + 3/2, 5/2; —). 


Both Wronskians can be essentially simplified by means of recurrency relations between 
Confluent Hypergeometric Functions: 


2 aF(a + 1, 3/2; = (2a - l)F(a, 3/2; ^f-) + F(a, 1/2; ^); 


3 E 2 


3E 2 


3 E' 2 


E F(a + 3/2, 5/2; —) = 4 (F(a + 3/2, 3/2; —) - F(a + 1/2, 3/2; —) ), 

and the following expressions are obtained: 

3 


W[& 1 \x,E),^ 1 \x)] +0 = 


| \ o rp2 o rp2 

—F(a, 1/2; ) + (2a - 1 )EF(a, 3/2; ) 


W[t' 2 >(x, E),^\x)] +0 = 3£(g 8 A) f(a + 1/2,3/2; 1|-) - F(a + 1/2,1/2; ^). 


After that, taking into account (1041) . the SGE (HOI) for U = —/, takes the form: 


(2a - 1 )yF(a, 3/2; y 2 ) + -(y + r/)F(a, 1/2; y 2 ) 


r(a 

—r(a + 1/2) 


y(y - y)F(a + 1/2, 3/2; y 2 ) - F(a + 1/2,1/2; y 2 ) 


= 0, 


where new variables are used: 


y = 


V3E 

2V2 ] 


n = 


V3\ 

2\/2 ’ 


(44) 


(45) 
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and the following constraint must hold: 


4 a = ( y 2 — y 2 ). (46) 

The obtained equation (l44]h from which one has to determine solutions y n , and therefore, 
by (j45l) - the allowed values of energy E n , is transcendental equation. In contrast to well- 
known exactly solvable models, these values can not be written as an explicit function of 
n. In this sense the model is similar to the standard quantum problem of a particle in 
the rectangular well of finite depth, which is considered traditionally in many textbooks. 
The only opportunity to solve such sort of problems is to use some appropriate numerical 
methods. 

In the present case, the suitable tool is to explore the computer algebra system MATH- 
EMATICA [27j. To simplify this task, it is technically useful to divide both sides of (1441) by 
a product r(a)r(o: + 1/2), this is equivalent transformation for a ^ —k/2 ; k = 0,1,2,..., 
where this product is infinite. One can fix some values for parameters to illustrate the solu¬ 
tion of SGE. The choice y — — 1 will not be considered here since it corresponds to parameter 
a = 1/2, leading to potential (1T81) without first term. Choosing, for example, y = —2 or 
y = —3, values of several lowest y n can be determined with good accuracy from the plot of 
l.h.s. of (|44|) after mentioned division by Gamma’s. The result is given at Figs.4-5, where 
the discrete values y n (proportional to energy eigenvalues E n ) are plotted as functions of the 
level number n for y = —2 and y = —3. 
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Fig.4 Discrete values y n as function of n for y = — 2 . 
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Fig.5 Discrete values y n as function of n for r/ = —3. 


One can notice that the first level n — 1 on both plots corresponds to y\ = r/, i.e. according 
to (J45]), to E\ — A. This observation is universal, since one can derive from the expression 
([3l]) that <F ( 1 ) (x) ~ T 0 (a;, a) for E = A and a = 0, and it has no nodes for finite x. For 
a = 0, relation (l34l) is fulfilled only if the constant A ^2 = 0. Comparing this argumentation 
with SGE, one can conclude that the second term in Eq. (l43l) vanishes, and therefore, Eq.(j44j) 
contains only first term. For a = 0 and y = ij this term is zero, and SGE is fulfilled. Thus 
the states with E\ — \ are the physical ground states in the case U = — I. Concerning first 
excited states, one can check by means of MATHEMATICA that for r) = —2 their wave 
functions with E 2l A 3 , K 4 from Fig.4 have just 1,2,3 nodes for finite x. All these levels for 
77 = —2, starting from the ground state E 2 can be seen at the Fig. 6 , where the l.h.s. of SGE 
(ITT) divided by product of Gamma functions is plotted. The positions of roots coincide with 
that on Fig.4. 



Fig.6 The plot of l.h.s. of (|44|) for rj = —2 divided by r(a)r(a + 1/2). Nodes correspond to 

the discrete spectrum y n . 
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The analogous procedure can be performed for the case U = +/. In this case, Wronskians 
with yd 2 \ instead of yd 1 ), must be calculated. The final result for the SGE is again the 
transcendental equation: 

w fofo + V 2 ; y 2 ) + 2 ^( 2a - l ) F ( a , 3 / 2 ; s/ 2 )]-^ 1 , 3 / 2 ; - ? / 2 ) - 

T(a + 1/2) [ 

—F(a, 1/2; ?/ 2 ) | - 

-Jr {Ml/ - ^(« + 1/2, 3/2; ?/ 2 ) - V F(a + 1/2,1/2; y 2 )]F(l, 3/2; - ?/ 2 ) + 

T(a) [ 

+yF(a + 1/2, 3/2; ?/ 2 ) j> = 0. (47) 

7. Conclusions. 


According to [18], [T9] , an arbitrary supercharge operator ([9]) allows factorization onto 
two first order supercharges with real or complex intermediate potential depending on the 
sign of constant d. For the pair of superpartner potentials (ITSlh factorization of second order 
supercharges ([9]) as in [18], [T9] gives the intermediate potential: 

T . . , f" 3 f' 2 A 2 A/' 2/3 7 

v int{x) -f- — + + — + — -X + + 


A 4a 
+ 


(48) 


36x 2 ' a; 4 / 3 3x 2 / 3 ’ 

which is real. It is obvious that no shape invariance between the systems with potentials 
and (g8D is observed. In this sense, the shape invariance between systems with potentials 
(USD is irreducible. The intermediate Hamiltonian with potential (ITS]) is intertwined with 
the initial H by means of first order operators of the form (jTJ) with superpotential 


W ( x , a) = f(x) - i- + 2 ^ x ^ ( 49 ) 

Necessary to mention here, that besides of [28], where mainly the reflectionless potentials with 
two-step shape invariance were considered, the shape invariance of second order was studied 
also in a few recent papers [29], [30]. The first of them deal with factorized supercharges and 
the superpotential W linearly depending on parameter a of shape invariance. In contrast to 
this approach, in the present paper in terms of W we considered actually the ’’square root 
dependence” of W via y— d = a fa (see m3). As for the [30], many potentials with two-step 
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shape invariance were obtained there but most of them can not be written explicitly - only 
in terms of rather complicate functions z(x). 

The main results of the present paper are the following. Using the shape invariance ap¬ 
proach in the framework of second order SUSY QM, a new one-dimensional shape invariant 
model was built, and its shape invariance is irreducible. The corresponding Schrodinger 
equation admits directly the exact solution. Due to the singular properties of the obtained 
potential, the model is well defined only after the suitable choice of the connection condi¬ 
tions at the singular point. These conditions were rewritten explicitly as specific Spectrum 
Generating Equations (SGE) for two simplest variants U = ±J. The SGE equations are tran¬ 
scendental equations for the discrete spectrum E n , which can be solved only numerically. 
Although MATHEMATICA is not able to calculate energies of very high excited states, the 
spectrum of the problem is evidently infinite in positive direction. The analysis of different 
terms in (|44|) in the context of their properties under reflection y —> —y shows that this 
equation having unbounded positive spectrum of y n can not have simultaneously very deep 
negative energies in the spectrum, confirming no fall to the center. 
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